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Abstract 

The NA60 experiment at the CERN SPS has studied low-mass muon pairs in 158A GeV In-In collisions. The mass 
and pt spectra associated with peripheral collisions can quantitatively be described by the known neutral meson 
decays. The high data quality has allowed to remeasure the electromagnetic transition form factors of the Dalitz 
decays //— >yu^yu"y and oj^p.^ pTif' . Using the usual pole approximation F - (\ - M^/A^)"' for the form factors, we 
find A-2 (in GeV^) to be 1.95±0.17(stat.)+0.05(syst.) for the 77 and 2.24+0.06(stat.)+0.02(syst.) for the co. While the 
values agree with previous results from the Lepton-G experiment, the errors are greatly improved, confirming now on 
the level of lOcr the strong enhancement of the w form factor beyond the expectation from vector meson dominance. 
An improved value of the branching ratio BR(a; p^^pTn^) - [1.73+0.25(stat.)+0.14(syst.)]-10"'' has been obtained 
as a byproduct. 
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1. Introduction 

The standard electromagnetic decay modes of light 
unflavored mesons (S=C=B=0) include the so-called 
Dalitz decays A —> Bl^T. Here, the meson A decays 
into an object B (a photon or another meson) and a lep- 
ton pair, formed by internal conversion of an interme- 
diate virtual photon with invariant mass M. Assum- 
ing point-like particles, the decay rate of this process 
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VS. M can exactly be described by QED fll]. How- 
ever, the rate is strongly modified by the dynamic elec- 
tromagnetic structure arising at the vertex of the transi- 
tion A — > B. This modification is formally described by 
a (multiplicative) transition form factor \FABiM)\^. A 
major element governing \Fab\^ is the resonance inter- 
action between photons and hadrons in the time-like re- 
gion, commonly referred to as vector meson dominance 
(VMD). Experimentally, \Fab{M)\^ is directly accessi- 
ble by comparing the measured invariant mass spectrum 
of the lepton pairs from Dalitz decays with the point-like 
QED prediction. A comprehensive review of the topic 
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is contained in iQl ■ 

The physics interest in studying DaHtz decays and 
the associated transition form factors is twofold. First, 
the electromagnetic interaction continues to be an ex- 
tremely useful tool to gain deeper insight into meson 
structure, while the role of the resonance interaction 
in this context is far from being quantitatively settled. 
Because the experiments are very difficult, the quality 
of the existing data is generally poor. Second, and re- 
lated to the last point, the study of direct production of 
dileptons in high-energy nuclear collisions in the con- 
text of thermal radiation requires a precise and com- 
plete knowledge of the characteristics and the relative 
weights for the existing decay channels, and this is uni- 
versally true at all facilities where such studies are on- 
going (SIS, SPS, RHIC, and FAIR in the future). Disre- 
garding the case of the if', the major two Dalitz decays 
contributing to the mass range M>0.2 GeV are those 
of the 77(548) and the w(782). For the dielectron chan- 
nel, the existing results on \F^ for the 77 and w (01 
are not accurate enough for meaningful physics conclu- 
sions. For the dimuon channel, however, significant re- 
sults on |Fp have been obtained by the Lepton-G exper- 
iment, both for 77 A^^/U T and for oj ji^ijrn'^ |0]. 
Using the usual pole approximation \2A 



\FY- = (1 -M7A0 



2\-2 



(1) 



for the form factors, A ^ has been found to be 1.9+0.4 
GeV-2 for the 77 and 2.36+0.21 GeV ^ for the u. While 
the value for the 77 is compatible with VMD within 
its large error, the value for the o) exceeds that ex- 
pected from VMD (1.69 GeV"^) by 3 standard devia- 
tions. This discrepancy, statistically significant, has re- 
mained unexplained up to today. Numerically, the asso- 
ciated enhancement of the mass-differential decay rate 
relative to that for VMD amounts to about one order 
of magnitude at M = 0.6 GeV, i.e. close to the kine- 
matic limit of - M„o - 0.648 GeV, with correspond- 
ing consequences for values and systematic errors of 
the yield of excess dileptons observed in this mass re- 
gion (H i i E [H - 

In this Letter, we present new results on the transi- 
tion form factors of the Dalitz decays 77 IJ^fJ^y and 
(1) — > ^^fTTfi. They have been obtained as a byproduct 
of the ongoing analysis of low-mass dimuon production 
in 158A GeV In-In collisions, exploiting here the nearly 
pp-\\ke peripheral rather than the more central interac- 
tions associated with excess dileptons iSllIil. The 
high data quality has enabled us to greatly improve the 
accuracy of the form factor measurements as compared 
to the Lepton-G experiment. 



2. Experiment 

The NA60 experiment at the CERN SPS is described 
in detail in 11211 . In short, the apparatus complements 
the muon spectrometer previously used by NA50 with a 
high-granularity radiation-hard silicon pixel telescope, 
placed inside a 2.5 T dipole magnet. The telescope 
tracks all charged particles upstream of the hadron ab- 
sorber and determines their momenta independently of 
the muon spectrometer. The matching of the muon 
tracks before and after the absorber, both in coordinate 
and momentum space, strongly improves the dimuon 
mass resolution in the low-mass region and reduces the 
combinatorial background due to n and K decays. The 
additional bend by the dipole field greatly improves 
the opposite-sign dimuon acceptance at low masses 
and low transverse momenta. The rapidity coverage is 
0.3<ycm<l-3 in this region. The selective dimuon trig- 
ger and the radiation-hard vertex tracker with its high 
read-out speed allow the experiment to run at very high 
rates for extended periods, leading to an unprecedented 
level of statistics for low-mass lepton pairs. 

3. Analysis Procedure 

The results reported in this Letter were obtained from 
the analysis of data taken in 2003 for 158A GeV In- 
In collisions. The analysis procedure is also described 
in detail in lll2ll . The essential steps of the data recon- 
struction concern the tracking in the two spectrometers, 
vertex finding, and matching of the tracks. Matching is 
done by selecting those associations between the muon- 
and pixel-spectrometer tracks which give the smallest 
weighted squared distance {matching x^) between the 
two tracks, in the space of angles and inverse momenta, 
taking into account their error matrix [12]. The com- 
binatorial background of uncorrected muon pairs orig- 
inating from n and K decays is determined by a mixed- 
event technique. After subtraction of the combinato- 
rial background, the remaining opposite-sign pairs still 
contain "signal" fake matches (associations of genuine 
muons to non-muon vertex tracks). These have a shape 
of the matching distributions different from those of 
the true matches. They are determined either by an over- 
lay Monte Carlo method (used here) or by event mix- 
ing [12], with identical results, and are then also sta- 
tistically subtracted from the data. The collision cen- 
trality of the events is defined through the total charged- 
particle rapidity density as measured by the silicon pixel 
telescope. 

For the purpose of this Letter, solely peripheral In- 
In collisions are considered. To keep sufficient events. 



they are selected through the cut in multiplicity density 
4<dNch/dTi<30, with an average multiplicity density 
{dNchldrj)-\l . The raw opposite-sign, background and 
signal dimuon mass spectra for this peripheral selection 
pre shown in Fig. [T] After subtracting the relatively 



while GEANT is used for transport through the accep- 
tance of the NA60 apparatus, including the effects of 
the dimuon trigger The Monte Carlo data are overlaid 
onto real data and then reconstructed in the same way 
as the latter, to take account of the pair reconstruction 
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Figure 1: Mass spectra of the opposite-sign dimuons (upper his- 
togram), combinatorial baclcground (dashed), signal fake matches 
(dashed-dotted), and resulting signal (histogram with en'or bars). 



small contributions from combinatorial background and 
signal fake matches, the resulting net spectrum contains 
about 26 000 muon pairs in the mass range 0.2- 1 .4 GeV. 
The average signal-to-background ratio is ~2.5. Al- 
though the relative uncertainties of the combinatorial 
background are ~4% for the peripheral selection (larger 
than the 1% achieved for more central collisions lll2ll '). 
the resulting systematic errors of the net data are still 
on the level of only about 1.5%. The vector mesons co 
and (j) are completely resolved; even the rare two-body 
decay rf — > ji^ii' is seen. The mass resolution of the lo 
is 20 MeV. 

As shown in our previous analysis flsl]. the peripheral 
data can fully be described by the expected electromag- 
netic decays of the neutral mesons. In the procedure 
used then and updated now, muon pair production from 
the 2-body decays of the //, p, u and (p resonances and 
the Dalitz decays of the rj, rj and co is simulated using 
the improved hadron decay generator GENESIS ll411 . 



Figure 2: Signal pairs after subtraction of the total background. A 
superposition of the known meson decays describes the data quantita- 
tively (see text). 

efficiency. The data are fit with this "decay cocktail" of 
sources, using the production cross section ratios rj/o), 
pi to, (pIcD and the level of dimuons from charm {DD) 
decays as freeparameters; the ratio rj ja) is kept fixed 
at 0.12 1I14I llSn . The branching ratios of the different 
decays are taken from the PDG \ and the transition 
form factors of the three Dalitz decays are those mea- 
sured by Lepton-G 10. 301 (which is also the source for 
the respective three branching ratios in the PDG). The 
Pt spectra of all hadrons {rj, u, (p and p) entering into the 
acceptance filtering have been precisely measured over 
the full /?r-range, including their centrality dependence, 
and those for the peripheral selection are used here. The 
y distributions are also measured, confirming those in 
the original code ifTill . The coie„ distributions {fics is 
the polar angle of the muon angular distribution in the 
Collins-Soper reference frame) are uniform for the reso- 
nance decays as recently measured fioli, and the angular 
distributions of the Dalitz decays are (l-i-cos'^) for the 



rj and uniform for the cj lll4ll (here 9 denotes the polar 
angle relative to the virtual photon direction). Data and 
fits, including an illustration of the individual sources, 
are shown in Fig. |2l The fit quality is good through- 
out. The only reminiscence to the excess dimuons found 
at the higher centralities fs', "9"] is a slightly enhanced 
pIlo ratio as compared to pp interactions 1 13], attributed 
to some contribution from n^n'' annihilation already at 
{dN ch I dT])-\l (without in-medium effects). The form 
factor anomaly of the cj Dalitz decay, responsible for 
the peculiar shape close to the kinematic cut-off in con- 
trast to the 77, is well visible in Fig. |2] and seems to be 
required by the data, although the description is not per- 
fect (we will come back to that later). It is important 
to add that the differential acceptance variations in the 
full M-pT plane, including a decrease by two orders of 
magnitude in the region of low M and low pj ifisll . have 
been understood to within <10% on the basis of the ob- 
served 777- -independence of the particle ratios extrapo- 
lated to full phase space jTsll . suggesting a significantly 
better accuracy in the mass domain alone. 

In the subsequent analysis, we will turn the proce- 
dure around. We will isolate the Dalitz decays of the 
rj and u as well as possible and measure the associated 
transition form factors without any a priori input to the 
description of the data in this region. The influence of 
all other decay sources on the results will be discussed 
in detail separately, considering the uncertainties con- 
nected to them as sources of systematic errors. 

In the first step, the 2-body decays of the narrow vec- 
tor mesons w and (p are subtracted in the same way as 
done before to isolate the excess dimuons at higher cen- 
tralities The yields are determined such as to 
get, after subtraction, a smooth underlying continuum. 
As discussed previously [j9|], the accuracy of this proce- 
dure is very high, about 3-4% for the cj and 2% for the 
0. These two sources are, in any case, completely out- 
side of the mass window relevant for the study of the 77 
and (JL> Dalitz decays, i.e. the window 0.2<M<0.65 GeV. 
The sole reason for their subtraction is the isolation of 
the broad vector meson p, which is normally masked 
by the much narrower u at nearly the same mass, mak- 
ing it then much easier to control the systematics due to 
a small contribution from the low-mass tail of the p in 
the mass region of interest here. Within that region, the 
well resolved 77 p^p' channel, with a mass resolution 
of about 13 MeV, is also subtracted with high accuracy, 
based on the same criterion as for the to and 4>- The 
remaining sources, the rj Dalitz decay and charm, are 
only on a level of a few % each of the total yield. As will 
be shown in next section, the final results are completely 
immune to the treatment of these sources, but they will 



ultimately also be taken out. After all subtractions, the 
remaining sample size is about 15 000 pairs, ~9 000 for 
the 77 Dalitz, ~3 000 for the w Dalitz and ~3 000 for the 
p. The corresponding Lepton-G numbers for the two 
DaUtz decays are 600 and 60, respectively 

The treatment of the acceptance of the NA60 appara- 
tus will also be turned around. Instead of dealing with 
the results after the acceptance filtering as in Figs.[T]and 
|2] we now correct the net data obtained from the first 
step for acceptance. The final physics outcome is, of 
course, invariant as to whether the analysis is done at the 
input or the output, i.e. before or after the acceptance fil- 
tering. The advantages of the reversal of the procedure 
are twofold. The large number of parameter variations 
in the fits to the residual data (>60) can much more ef- 
ficiently be done at the input, without the need to prop- 
agate each choice through the complete Monte Carlo 
chain. In addition, the final results can then be judged 
on the basis of the original spectral shapes, without ac- 
ceptance distortions. The acceptance is determined for 
the mixture of the three sources left in Fig.|3] i.e. the 77 
and (jj Dalitz decays plus the p, using precisely the same 
input distributions as before and the weights as obtained 
from the fits to the peripheral data (Fig. |2]i. The result- 
ing mass dependence is rather flat, with variations in 
the region of the Dalitz decays, 0.2<M<0.65 GeV, by 
only a factor of 2. A rise by a further factor of 2 occurs 
across the mass region of the p. For systematic studies, 
the relative weight of the iw-Dalitz decay in the mixture 
is varied by a factor of 1.8 (for reasons to be discussed 
later), resulting in local changes of the acceptance by 
only <+3%. 

4. Results 

The net mass spectrum of the muon pairs after sub- 
traction of the three narrow resonances 77, o) and ^ as 
well as the 77 Dalitz decay and charm (see below), cor- 
rected for acceptance and pair efficiency, is shown in 
Fig. |3] The ordinate is in a.u. and does not any longer 
reflect the measured number of counts. The spectral 
shape of the data in this figure looks impressive. Be- 
yond the 77 Dalitz decay which was easily recognizable 
before in Figs.[Tl|2]due to its dominance in the mass re- 
gion M<0.5 GeV and its characteristic mass shape, the 
p now shines out completely isolated, and the o) Dahtz 
decay in between becomes directly recognizable in the 
mass window 0.5<M<0.65 GeV through the character- 
istic shoulder close to the kinematic cut-off of its mass 
distribution. This shoulder reflects, beyond any doubt, 
the qualitative existence of the strong anomaly in the 
associated electromagnetic transition form factor QSl^ 



and the data quality in this well-isolated section raises 
the expectation that the extraction of the quantitative de- 
tails will be possible with a high reliability. This is in- 
deed the case. 
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Figure 3: Acceptance-corrected mass spectrum of the muon pairs after 
subtraction of the oj, resonance decays and the nearly neghgible 
contiibution of the Dahtz decay and charm. The hnes represent a 
fit with Eq. (sj, showing the individual contributions from the rj and 
OJ Dalitz decays, the p resonance decay and the sum of the three (see 
text for final fit parameters used here). The shaded areas indicate the 
expectations for the f; and a) Dalitz decays for point-like particles, 
defined by [QED] 1 1] in Eqs. <2l3t . 



The further study is based on global as well as lo- 
cal fits to the data, using a superposition of the three 
sources with (mostly) 5 free parameters: the yields of 
the three sources and the pole parameters of the form 
factors of the two Dalitz decays, defined already above 
with Eq. ([T]). All other physics parameters which are not 
strictly known will be varied to study their influence on 
the fits. Once the final fit parameters and their errors are 
fixed, the three sources in Fig. [3] will also be disentan- 
gled, making it then possible to present the form factors 
|F,(M)p of the two Dalitz decays in the usual way. 

Analytically, the dilepton mass spectrum of the Dalitz 



decay // fi^fi y is given by yHSfl (M = m^^) 
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while that for the Dalitz decay w — » yU^// ifi is defined 
by the slightly more compUcated expression 10,01 
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For the resonance decay p — » n^iT, we use the line 
shape (including a Boltzmann term) characteristic for 
hadro-production of the p ifTiiflTll 
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The fit function with the 5 free parameters A, and A. 
can then be written in the form 
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■ AJ,{M, A,:2) + A^UM, + Apfp(M) (5) 



where the mass differential functions f, contain the full 
information of Eqs. (|2|-|4]l. 

The specific fit shown in Fig. |3] illustrates already the 
final outcome: the individual fit lines for the two Dalitz 
decays are above the respective QED expectation, im- 
plying the form factors to be consistently > 1 and actu- 
ally » 1 for the (L) spectrum close to its kinematic limit. 

The interest in leaving the yields of the Dalitz de- 
cays also free in the fits is connected to their branch- 
ing ratios. As a byproduct to the form factor studies, 
these can directly be measured by relating the yields to 
those of the respective resonance decays, which are sub- 
tracted but known. In case of the rj, the experimental 
error of the decay i] —> fi'^jj." is unfortunately ~30%, 
much larger than those of either of the two branching 
ratios [16]; there is thus no chance for improvements. 
In case of the co, however, the errors of the reference 



decay co — > n^fJT are quite small: 3-4% for the mea- 
sured yield in Fig. |2] and 2% for the branching ratio, if 
we replace the ji^^^ value (9.0+3.1)- 10"^ by the much 
more accurate e^e" value (7. 16+0. 12)- 10"'' |16] as jus- 
tified by lepton universality. In contrast, the error of the 
branching ratio of the Dalitz decay BR(w — > yu^yu tt*') 
= (9.6+2.3)- 10"^ [6, 16] is much larger, 24%, opening 
a realistic chance for a more accurate measurement of 
the latter The outcome will indeed be a larger value, 
and the two acceptance options mentioned before take 
care of this initial ambiguity in the mixture of sources, 
labeled in the following A (for the new value) and B (for 
the PDG value). 

Altogether, about 60 different fits were done, varying 
a number of parameters and their combinations. The 
parameters were the following. The acceptance had the 
options A and B. The temperature parameter Tp of the 
p, so far unmeasured in /7/7-like interactions, had the op- 
tions 170 and 140 MeV, while the pole mass and width 
of the p were (mostly) fixed at Mp - Q.llQ and Tp - 
0.150 GeV flu], respectively. Since the acceptance- 
corrected data still contain the resolution smearing of 
the NA60 set-up, resolution smearing was imposed on 
the fit function Eq. (|5]l to study the sensitivity. The con- 
tribution of the Dalitz decay was either fully left in the 
data sample, or subtracted under the different assump- 
tions tj'/lo^O.U 01 m, 0.24, 0.36 or 0.48, amount- 
ing to subtracted fractions of only 1%, 2%, 3% and 
4% of the full yield in the mass region M<0.64 GeV. 
The contribution from charm {DD decays) was either 
fully left in the data sample, or subtracted under the dif- 
ferent assumptions of 30%, 60% or 100% of the total 
yield in the mass window 1.2<M<1.4 GeV to be charm 
(and Drell-Yan), amounting to subtracted fractions of 
1%, 2% and 4% of the total yield in the mass region 
M<0.64 GeV. Specific fits were done (i) with a branch- 
ing ratio of the oj Dalitz decay frozen to the PDG value, 
and (ii) with the temperature parameter Tp of the p also 
left as a free parameter Fit ranges were mostly global, 
covering the complete mass range 0.2<M<0.9 GeV, but 
for specific goals also local, like 0.5<M<0.9 GeV. The 
quality of the fits for each parameter combination was 
judged by the respective x^/ndf, assessed globally in 
the complete mass range 0.2<M<0.9 GeV, but also lo- 
cally in the subwindows 0.2<M<0.48 GeV (-80% t] 
DaUtz), 0.48<M<0.66 (-67% oj DaHtz), 0.66<M<0.86 
(-100% p), and 0.75<M<0.95 GeV (upper tail of the 
p). The detailed results from the individual 60 fits, in- 
cluding the values of the fit quantities, their errors, the 
X^/ndf in the individual windows and a number of fur- 
ther aspects of the analysis can be found in a publicly 
accessible NA60 Internal Note I.18J . Here, we restrict 



ourselves to a summary of the results and the associated 
systematics. 

Very generally, the fits have been found to give re- 
markably stable and reproducible results. The fit quality 
in the mass region of the two Dalitz decays, i.e. M<0.65 
GeV, is completely insensitive to variations of the ac- 
ceptance, of Tp, of the resolution folding, and of the 
fraction of subtraction of the 77 and charm contributions. 
The values of x^/ndf are always -1, and the variations 
of the extracted fit parameters are mostly <l/2 of the 
(statistical) fit errors. Given this situation, we quote the 
final values of the fit parameters as the (unweighted) av- 
erage over all measured values, their statistical errors as 
the average over the fit errors (which hardly vary at all), 
and their systematic errors as the rms deviation of the 
individual values from the average. The ratio system- 
atic/statistical errors is about 0.3 for the two pole param- 
eters, and about 0.5 for the w-Dalitz branching ratio. 

However, the mass region outside of interest here, i.e. 
M>0.65 GeV and in particular the high-mass tail of the 
p, is indeed sensitive to the variations of some of the pa- 
rameters, and the extracted values of x^/ndf can reach 
up to values of 7. The global conclusions from this part 
of the analysis are a clear preference for the higher tem- 
perature of the p, for a subtraction of the 77 Dalitz con- 
tribution on the level of at most rj /w=0.12 iHlil, and 
for a full subtraction of the charm contribution. It is for 
this reason and not for reasons of any sensitivity in the 
region of the 77 and w Dalitz decays, that these specific 
77 and charm contributions have been subtracted for the 
data sample selected for Fig. [3] 

A final comment on systematics. The analysis proce- 
dure as used is self-consistent; the acceptance has been 
assessed on the basis of measured data throughout, in- 
cluding the anomalous form factor of the oj. However, 
the results are extremely robust as to deviations from 
that. For the a>, e.g., usage of the VMD form factor in 
the simulations leads to drastic deficits in the descrip- 
tion of the (acceptance-filtered) data in Fig.|2]fl8ll. Yet, 
an acceptance correction based on such an inferior de- 
scription, with quite different (mass-dependent) weights 
of the 3 sources to before, still leaves the characteris- 
tic shoulder of the a> Dalitz decay (Fig. O essentially 
unchanged, and the fit value of the pole parameter is 
found to only change by about 1 standard deviation (sta- 
tistical) [18]. It hardly needs to be stressed that the 
self-consistent results obtained from fits at the input, i.e. 
after acceptance correction, are found to be absolutely 
identical to those obtained from fits to the directly mea- 
sured data, i.e. before acceptance correction il8ll . 

In detail, the following numerical results have been 
obtained. 



The pole parameter of the electromagnetic transi- 
tion form factor of the Dalitz decay rj — > ju^// y 
is measured to be A^2=1.95+0.17(stat.)+0.05(syst.) 
GeV"^. It perfectly agrees with the previous mea- 
surement of the Lepton-G experiment A,^^= 1.90+0.40 
GeV"^ as well as with predictions from VMD, A,^^=1.8 
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Figure 4: Experimental data on tlie ;;-meson electromagnetic transi- 
tion form factor (red triangles), compared to the previous measure- 
ment by the Lepton-G experiment (open circles) and to the expecta- 
tion from VMD (blue dashed Une). The solid red and black dashed- 
dotted Hnes are results of fitting the experimental data with the pole 
dependence Eq. Q. The normalization is such that \F^(M = 0)|=1. 

GeV"^ The characteristic mass A is equal to 

A^=0.716+0.031(stat.)+0.009(syst.) GeV, as compared 
to the value from Lepton-G of A,,=0.724+0.076 GeV 
or to the VMD value of A;,=0.745 GeV. Our result im- 
proves the Lepton-G error by a factor of 2.3, equivalent 
to a factor of 5 larger statistics. The error improvement 
to be expected from the difference in sample sizes (9 000 
vs. 600) would have been larger (a factor of 3.8), but 
this is only found if the cj Dalitz decay is frozen in the 
fitlH. 

The pole parameter of the electromagnetic transi- 
tion form factor of the Dalitz decay lj ju^fi'n^ 
is measured to be A^^ - 2.24+0.06(stat.)+0.02(syst.) 
GeV"^. Within errors, it agrees with the Lepton- 
G value of A^2=2.36+0.21 GeV"^. Both experimen- 



tal results differ from the expectation of VMD of 
A„^=L68 GeV"^ [2]. The anomaly is therefore fully 
confirmed. The characteristic mass A is found to be 
A^=0.668+0.009(stat.)+0.003(syst.) GeV, as compared 
to the value from Lepton-G of A„ =0.65 +0.03 GeV or to 
the VMD value of A„=Mp=0.770 GeV. The confirma- 
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Figure 5: Experimental data on the w-meson electromagnetic transi- 
tion form factor (red triangles), compared to the previous measure- 
ment by the Lepton-G experiment (open circles) and to the expecta- 
tion from VMD (blue dashed line). The solid red and black dashed- 
dotted lines are results of fitting the experimental data with the pole 
dependence Eq. (T). The normalization is such that \F^{M = 0)|=1. 

tion of the anomaly receives particular weight through 
the fact that the statistical errors are improved by a fac- 
tor of nearly 4, equivalent to a statistics larger by a factor 
of >10. Referred to A"^, the previous measurement dif- 
fered by three standard deviations (3cr) from the VMD 
expectation, while our new measurement differs by lOcr. 
The error improvement to be expected from the differ- 
ence in sample sizes (3 000 vs. 60) would have been 
still larger (by a factor of 7), but thjs is only found if the 
7] Dalitz decay is frozen in the fit [18]. 

The branching ratio of the o) Dalitz decay BR(w — » 
jx^^'TiP) is found to be larger by a factor of 
1.79±0.26(stat.)±0.15(syst.) than that of the PDG \\^, 
i.e. Lepton-G [6], corresponding to a new absolute 
value of (1.72+0.25(stat.)+0.14(syst.))-10"^ Taking ac- 



count also of the 24% error of the PDG value, the two 
values of the branching ratio differ by 2 standard devi- 
ations, hardly significant. It is interesting to note that 
there is no contradiction between the increased branch- 
ing ratio as obtained from the present systematic fits 
and the previous description of the data on the basis of 
the PDG value as shown in Fig. |2] On the contrary, a 
deficit of the sum of the generated events in the mass 
region 0.48<M<0.64 GeV compared to the data (sub- 
tracting the J] — > ju^ju" channel also there) clearly ex- 
ists, reflected by ax^/ndf of about 5 in that region of 
Fig. |2] For M>0.55 GeV, this is partially compensated 
for by the effects of the higher pole parameter value 
from Lepton-G. Specific fits done to scrutinize such a 
compensation fail for the present data, however, thanks 
to the much higher data quality. Fixing the branching ra- 
tio to the PDG value and varying the pole parameter, the 
description of the data is already visually unacceptable 
for any choice of the pole parameter, and the minimal 
value of I ndf reached is about 3.3 1 1 8 ] . 

From a series of specific fits, leaving the tempera- 
ture parameter Tp in the line-shape description of the 
p also free, the average value of 7), is found to be 
170+19(stat.)+3(syst.) MeV. Since the event selec- 
tion is peripheral In-In, with {dNch/dT]}~ll , there will 
hardly be an in-medium influence on the value. We in- 
terpret this number to reflect the eff'ective temperature of 
the system at the time of creation, making it consistent 
with the same value of 170 MeV obtained by statisti- 
cal model fits of particle ratios in pp interactions. This 
is the first time in the literature that Tp has been deter- 
mined experimentally. 

The specific simulation lines shown in Fig. [3] are ac- 
tually based on the average fit parameters discussed in 
the preceding four paragraphs. Starting from here, we 
finally present the results on the transition form factors 
in the traditional way. In a first step, we isolate the in- 
dividual Dalitz contributions in the spectrum of Fig. |3] 
subtracting the contribution of the p ^^fT decay and 
disentangling the tj i-i^i^^j and co — > fx^ju^jfi decays 
as determined by the fits. This implies to use the same 
individual data points for the 77 and the co, subtracting 
for the rj the fit results of the <y and vice versa. With 
|F,(M)p— > 1 for M— >0, the individual normalizations are 
automatically fixed, and the QED and the form factor 
parts can be separately assessed. In a second step, the 
squared form factors |F,(M)p are obtained by dividing 
the difference data for the respective decay by its QED 
part. 

The results for the 77 and the lj are plotted in Figs. |4] 
and |5] respectively, keeping the data-point errors from 
Fig- El The pole parameters and theu- errors as obtained 



from the combined fits to both Dalitz decays are shown 
as inserts. Note that these are the correct values, while 
independent fits through the data points of Figs. |4] and 
|5] would automatically result in somewhat smaller er- 
rors of the pole parameters, since the respective other 
decay appears as fixed. The errors shown are statistical 
errors. The systematic errors are smaller than the statis- 
tical ones as outlined above. For the lj form factor in 
the mass region <0.45 GeV, where the 77 Dalitz decay 
dominates in the total mass spectrum (see Fig.O, this 
has explicitly been verified (on top of all other sources) 
by varying the form factors of the two Dalitz decays 
in the global fit procedure preceding the isolation JlSIl. 
Both figures also include the Lepton-G data 10, \Km 
and the expectations from VMD |0 El 01 for compari- 
son. Within the large errors of the Lepton-G data, per- 
fect agreement between the two data sets is seen in both 
cases, while the great improvement in data quality of 
the present results is completely apparent. Irrespective 
of the much reduced errors, the form factor of the 77 is 
still close to the expectations from VMD. The form fac- 
tor of the o), on the other hand, strongly deviates from 
VMD, showing a further (relative) increase close to the 
kinematic cut-off by a factor of ~10, and a factor of al- 
together ~100 relative to the pure QED part. 

Theoretically, the most elementary description of the 
transition form factors in terms of VMD does not only 
work reasonably well for the Dalitz decay of the 77, but 
also for that of the 77 , at least within the very large er- 
rors there |2]. The anomaly of the lj case has been a 
puzzle from the beginning. Slight enhancements of the 
oj form factor beyond VMD have been obtained histori- 
cally on the basis of a modified p propagator [T?] and a 
nonlocal quark model lEoll . More recent calculations on 
the basis of an effective Lagrangian approach to vector 
mesons ['21'], and of an extended VMD model including 
up to two excited p states ll22ll . do somewhat better, but 
now overestimate the form factor at low M and still un- 
derestimate it at high M in the region of the kinematic 
cut-off. One might rightly expect that our improved data 
will initiate new theoretical efforts to finally understand 
the physics behind. 



5. Conclusions 

To summarize, we have been able to measure the 
electromagnetic transition form factors of the 77 and at 
Dalitz decays with a much better precision than reached 
before, confirming after nearly 30 years the strong 
anomaly associated with the oj. A satisfactory theoreti- 
cal understanding is still pending. On purely empirical 



grounds, the new results will greatly diminish the uncer- 
tainties of direct dilepton measurements in this particu- 
lar mass region. As a byproduct, we have also obtained 
an improved value for the branching ratio of the lj Dalitz 
decay. 
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